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ABSTRACT 

We present the A-band local luminosity function derived froma sample of 
1056 bright (A < 15) A-selected galaxies from the Hawaii-AAO A-band redshift 
survey. The Hawaii-AAO A-band redshift survey covers 4 equatorial fields with 
a total area of 8.22 deg 2 . We derive both the non-parametric and Schechter 
luminosity function from our data, and determine M*(K) = —23.70 ± 0.08 + 
51og 10 (/i), a = -1.37±0.10 and 0* = 0.013 ±0.003 h 3 Mpc~ 3 for a universe with 
Q m = 0.3 and = 0.7. We also measure the A-band luminosity function for the 
early- and later-type galaxies from our morphologically classified subsample. It 
appears that later-type galaxies have a fainter M* and a steep slope, while early- 
type galaxies have a much brighter M* and a quite flat slope in their A-band 
luminosity functions. This is consistent with what have been found in optical 
type dependent luminosity function. The A-band luminosity density derived 
using our luminosity function is now measured at a similar redshift depth to 
optical luminosity densities in the SDSS redshift survey. It is 2 times higher 
than the previous measurement from the shallower 2MASS sample and resolves 
the previously reported discrepancies between optical and near-IR luminosity 
densities. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
Survey — galaxies: Near Infrared 
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1. Introduction 

The galaxy luminosity function is an important quantity in the study of galaxy evolu- 
tion and formation. Traditionally, galaxy luminosity functions have been derived in optical 
bands. It is now clear that extragalactic studies in optical bands suffer several systematic 
uncertainties and complexities compared to those using near infrared bands. In particular, 
dust extinction has little effect on K-band magnitudes; the K-correction in the i^-band is 
a much smaller and better understood quantity than in the optical bands, and it is inde- 
pendent of galaxy spectral types for z < 1. Because of this it is easier to detect a high 
redshift elliptical galaxy in an infrared band than in optical bands. A galaxy's near infrared 
luminosity is also a good tracer of its stellar mass independent of spectral type (Cole et 
al. 2001; Glazebrook et al. 2002). Recent theoretical studies show that the X-band galaxy 
luminosity function is a powerful constraint on galaxy formation theory (Baugh et al. 1998; 
Kauffmann k Chariot 1998). 

Most near infrared surveys, however, have been modest in size due to the small size of 
available infrared detectors. The advent of large format infrared array detectors has made 
possible a variety of wide-field near-infrared surveys, ranging from several at the 10 deg 2 - 
level (Gardner et al. 1996; Huang et al. 1997) up to the largest of them all, the 2 Micron All 
Sky Survey (2MASS) (Skrutskie et al. 1997). Obtaining optical redshifts for a K-selected 
sample, however, is difficult because the wide range in optical-infrared colours results in the 
requirement for a wide range of exposure times to acquire optical spectroscopy. In particular, 
very long exposures are required to secure redshifts for the reddest objects. Because of this 
there are only a few X-band luminosity functions available. The early X-band luminosity 
functions were derived from small size samples with number of galaxies ranging from 100 
to 500 (Mobasher et al. 1993; Glazebrook et al. 1995; Gardner et al. 1997; Szokoly et al. 
1998; Loveday 2000). After the second incremental release of 2MASS data, two teams (Cole 
et al. 2001; Kochanek et al. 2001) used the overlap between the 2MASS Extended Source 
Catalog (http://pegasus.phast.umass.edu) and two existing optical redshift databases, CfA2 
(De Lapparent et al. 1988) and 2dFGRS (http://www.mso.anu.edu.au/2dFGRS), to obtain 
very large (> 4000), if very shallow (K < 13), X-selected redshift samples to derive the local 
K-band luminosity function. However, the mean redshift for these samples is shallow, 0.025 
for the CfA sample(Kochanek et al. 2001) and 0.05 for the 2df sample(Cole et al. 2001). 

The observational goal of the Hawaii- AAO K-band redshift survey was to obtain a large 
medium-deep X-selected galaxy sample with redshifts and optical-infrared colors. Both 
optical (B and /) and near-infrared (K) images were taken at Mauna Kea Observatory 
with total infrared coverage of totaling 8.22 deg 2 (Huang et al. 1997). The spectroscopic 
observations were carried out on the Anglo- Australian Telescope (AAT) with the Two Degree 
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Field facility (2dF). In this paper we report K-band luminosity functions derived from a sub- 
sample of 1056 bright (K < 15) galaxies. The median redshift of this sample is 0.138 with 
the redshift distribution tail extending to z=0.5. In this paper, we adopt the Q m = 0.3 
and £7a — 0.7 and H = 100 km s -1 Mpc~ 3 cosmology model. A comparison of the if-band 
luminosity functions derived using different cosmology models is also presented. We briefly 
summarize the observation and data reduction in §2; in §3, we present the K-band luminosity 
functions; in §4, we compare our luminosity function with other's and discuss the difference; 
we summarize our results in 85. 



2. Observations 

Huang et al. (1997) have described in detail the acquisition and processing of our B-,I- 
and K-band imaging. In brief, the imaging survey was carried out at Mauna Kea Ob- 
servatory and represents over 70 nights of time on the University of Hawaii 88-inch and 
24-inch telescopes. The survey consists of four equatorial fields covering a total area of 
8.22deg 2 . The optical and near-infrared limiting magnitudes are Bn m it = 22, lu m it = 20.5, 
and Kumit = 16.0. The spectroscopic observations for the ^-selected sample were carried 
out at the AAT with 2dF, a multi-fiber spectrograph that can observe 400 targets within a 2 
degree field of view. These observations took approximately 10 clear nights spread out over 
the period 1997-1999. The exposure time for a spectroscopic observation was determined 
by the / magnitude of each object, for objects with / < 18.5 (I — K < 2.5 the K = 16 limit) 
1 hour exposure times sufficed, for the redder objects (18.5 < I < 21) exposure times of up 
to 4 hours were used. Weather variations meant that some 2dF observations were not as 
deep as others. When a redshift was not secured at a first attempt the object was flagged 
for re-observation in a later 2dF configuration. Multiple 2dF configurations were observed 
on each field to ensure maximal completeness. 

Our K < 15 sample is highly complete in redshift. From a total of 1201 objects imaged 
with K < 15, redshifts were obtained for 1056 galaxies. A further 59 objects turned out to 
be stars, leaving only 86 objects without identification. These objects are assumed to be 
galaxies for the purposes of incompleteness correction. Fig. 1 shows the redshift distribution 
of the sample. The median redshift for the K < 15 sample is 0.136. We calculate absolute 
magnitudes with a Q m = 0.3 and Q\ — 0.7 cosmology and use recently published K-correction 
due to Mannucci et al. (2001). K-corrections have also be tested using spectral evolution 
models which best fit the observed colours (Fioc & Rocca-Volmerange 1997). No significant 
difference was found in either the K-corrections or the resulting luminosity functions. 
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3. The Luminosity Function 



We use the Stepwise Maximum Likelihood (SWML) (Efstathiou et al. 1988) to derive the 
non-parametric luminosity function, and use the STY method (Sandage, Tammann & Yahil 
1979) to fit the Schechter function to the data. After obtaining M* and a by fitting Schechter 
function to our data, we use the minimum variance estimator (Davis & Huchra 1982) to 
determine the </>*. Both STY and SWML methods are based on the maximum likelihood 
principle, and therefore only valid for a complete galaxy redshift sample. Correction for 
our sample's (small) incompleteness is therefore needed. Zucca et al. (1994); Lin et al. 
(1996, 1999) have shown that a weighted STY method can be used to derive a luminosity 
function from an incomplete sample. In the weighted STY method, if a galaxy in the sample 
has no redshift, similar galaxies with redshifts (ie. similar magnitude, or similar colour, or 
both) receive increased weight in calculating the likelihood to represent the galaxy without 
a redshift, 



Here the p is the probability for galaxy i at redshift Zi, ip(L) is the luminosity function, 
and Wi is the weight for galaxy i. Wi = 1 for a complete sample, otherwise Wi > 1. 

In our case only 7% of targets do not have redshifts, so such a correction does not 
make a large difference. We fit Schechter functions to the sample both with, and without, 
correction. We follow Lin et al. (1996, 1999) in estimating the weight function for each galaxy 
with redshift in the sample. We divide the sample into apparent magnitude and color I — K 
bins, and the weight assigned to galaxy i with redshift is the ratio of total number of 
galaxies over the total number galaxies with redshifts in the bin where the galaxies i is. Since 
the incompleteness of the 4 fields varies, we have to calculate Wi for galaxies in each field 
separately. Table 1 summarizes the result of the STY fitting with 3 different cosmological 
models, and Fig. 2 shows the non-parametric K-band luminosity function and the Schechter 
function in the Vt m = 0.3 and Q\ = 0.7 universe. 

As shown in Table 1, the parameters of the Schechter function depend on the adopted 
cosmological models. This is because the median redshift is 0.138, where the effect of the 
geometry for the adopted cosmology model appears in calculating distance, and the change 
of a is due to the correlation between M * and a. There is a slight difference between the 
luminosity functions derived with correction and without correction. The M* calculated with 
correction is slightly brighter than that calculated without correction. This is because most 
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of the galaxies without redshifts are near the K — 15 limit and are either bright galaxies at 
relatively high redshifts comparing with rest of the sample, or local lower luminosity galaxies. 
Similar galaxies with redshifts in the sample are assigned more weight. The consequence is 
to cause M* to become somewhat brighter and the faint end stope slightly steeper after 
correction. 

Using the higher spatial resolution /-band images, Huang et al. (1998) was able to make 
a morphological classification for the very bright galaxies (K < 14). Among 225 galaxies 
with K < 14 in the sample, Huang et al. (1998) classified 111 galaxies as E/SO. Most of the 
remaining 114 galaxies are spirals, and only a few with peculiar morphologies are possible 
mergers or irregular galaxies. Therefore, we are able to derive the luminosity function of early 
type (E/SO) and later type (Spiral and Irregular galaxies) for the K < 14 sample. Fig 3 
shows the the K-band luminosity function for both types with the best fitting Schechter 
functions. The K-band luminosity function for the early type galaxies has a bright M* of 
M* = —23.56 ±0.26 and a flat slope a = —1.04 ±0.31, while that for the later type galaxies 
has a much fainter M* of M* = -23.28 ± 0.28 and a steeper slope a = -1.42 ± 0.31. 



4. Discussion 

Recently, several K-band luminosity functions were derived either directly from K- 
band surveys (Glazebrook et al. 1995; Gardner et al. 1997; Szokoly et al. 1998), or from 
optical samples with K-band images (Mobasher et al. 1993; Loveday 2000). These luminosity 
functions are roughly in agreement with each other, given that most samples are small in 
size and have large statistical uncertainty. To compare the K-band luminosity functions 
in different cosmology models, Cole et al. (2001) convert all existing i^-band luminosity 
functions to those in the Q m = 0.3 and = 0.7 cosmology model. Table 2 shows the 
parameters of Schechter function for all K-band luminosity functions, including those we 
adopted from Cole et al. (2001) and ours. We also list sample size and limiting magnitude 
for each survey. 

By comparing the Schechter functions, our luminosity function has a brighter M* and a 
steeper slope, only those of Szokoly et al. (1998) are close to ours. The rest of the luminosity 
functions have a M* at least 0.3 magnitude fainter than ours and a flat slope (a ~ — 1), 
including two 2MASS luminosity functions (Cole et al. 2001; Kochanek et al. 2001) and that 
of Gardner et al. (1997) who had a similar sample to ours. For more accurate comparison, we 
plot the luminosity functions against each other in Fig. 4. Fig. 4 shows that our luminosity 
function is higher and steeper than the 2MASS luminosity functions and that of (Gardner et 
al. 1997). This is in agreement with the comparison in Schechter functions that the 2MASS 
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luminosity functions have fainter M*, lower (ft*, and flat a. 

There are two different important aspects between our sample and the 2MASS sam- 
ple(Cole et al. 200f ; Kochanek et al. 2001): different ways of measuring total magnitudes for 
galaxies and different redshift ranges for both sample. Each of them could cause the differ- 
ence in deriving luminosity function. If we assume that both photometric methods measure 
true total magnitudes for X-selected galaxies, the difference between both i^-band luminos- 
ity functions implies a change of luminosity function at different redshift ranges. Since our 
sample covers a wider redshift range, we are able to derive the i^-band luminosity function 
in the lower redshift bin. This will allow us to test if we are able to reproduce the results 
from the 2MASS sample using our lower redshift bin sample. The limiting magnitude in 
both Cole et al. (2001); Kochanek et al. (2001) samples are around K ~ 13, which is too 
bright for our sample to be statistically significant. However we can adopt a redshift limited 
subsample which approximates their redshift range, it will simply extend ~ 2 magnitudes 
further down the luminosity function. Thus we set the limiting redshift z < 0.1, to approx- 
imate the redshift range to those of the 2MASS sample while allowing us to have enough 
galaxies to derive the luminosity function. Fig. 4 shows that our luminosity function for the 
subsample with z < 0.1 matchs those luminosity functions of (Cole et al. 2001; Kochanek et 
al. 2001) very well. We also fit the Schechter function to our subsample with z < 0.1, and 
obtain that M* = -23.10 ± 0.15, a = -0.93 ± 0.16, and (ft* = 0.012 ± 0.004 Mpc~ 3 consis- 
tent with those derived using the 2MASS samples. This implies that the normlisation of the 
luminosity function is a function of redshift, even at lower redshifts, and that this explains 
the apparent discrepancy in published estimates without needing to invoke systematic errors 
in deriving total magnitudes in both samples. 

The steep slope of our K-band luminosity function is actually consistent with those of 
optical luminosity functions obtained in several current wide-field optical surveys. Folkes et 
al. (1999) drived a 5-band luminosity function with a = —1.28 ± 0.05 using a large redshift 
sample obtained in the 2df Galaxy Redshift Survey (2dFGRS). Blanton et al. (2001), and 
the SDSS team obtained the optical luminosity functions in u*, g*, r*, i*, and z* bands, and 
the a for these luminosity functions in the 5 bands are —1.38, —1.26, —1.20, —1.25, and 
— 1.24 respectively. Cole et al. (2001) used the average color to transfer the SDSS z*-band 
luminosity function to a K-band luminosity function with M* = —23.67, (ft* = 0.0127, and 
a = —1.24, which are very close to the parameters of our luminosity function. We would like 
to point out that the peak of the redshift distribution for the SDSS sample is about 0.1, much 
closer to ours than to those of the 2MASS sample (Cole et al. 2001; Kochanek et al. 2001). 
The implies that the normalization of the luminosity function is a function of redshift, even 
at lower redshifts, and that this explains the apparent discrepancy in published estimates 
without needing to invoke systematic errors in deriving total magniutdes in both sample. 
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Before the 2MASS sample was available, there were no K-selected samples with mor- 
phological classification, hence no type-dependent K-band luminosity functions. Loveday 
(2000) derived K-band luminosity function for the emission line galaxies (ELG) and also for 
the galaxies without emission lines (Non-ELG). He found that the M* of the K-band Lumi- 
nosity function for ELG is one magnitude fainter than that for non-ELG. For the first time, 
Kochanek et al. (2001) was able to derive the i^-band luminosity function for both early 
and later type galaxies. We can compare ours with those of Kochanek et al. (2001). Since 
both samples have very bright limiting magnitudes, K < 11.2 for Kochanek et al. (2001) 
and K < 14 for our sample, the effect of the geometry for the adopted cosmology models 
cannot make any significant difference in calculating absolute magnitudes at such a low red- 
shift. We can compare the two results directly. In Table 3 we list the M* and a derived 
from both our K < 14 sample and the 2MASS sample (Kochanek et al. 2001). Our if-band 
luminosity function for early type galaxies has a similar M* and a to those of Kochanek et 
al. (2001), but the luminosity functions for later type galaxies are more different: our M* is 
0.3 magnitude brighter and our slope is much steeper (a = —1.42) than theirs (a = —0.87). 
Steep slopes (a < —1) for later type galaxy luminosity functions are also seen in the optical 
bands (Marzke et al. 1994; Bromley et al. 1998; Marzke et al. 1998; Folkes et al. 1999). 

The spectrum of the universal luminosity density is another way of testing consistency 
between optical and infrared luminosity functions (Wright 2001). Dwek et al. (1998) sug- 
gested that the spectrum of the universal luminosity density can be fitted by an average 
spiral galaxy SED (Schmitt et al. 1997). Wright (2001) used the average spiral galaxy SED 
model to fit the luminosity densities derived from the SDSS luminosity function, and find 
that the model predicts the luminosity densities at J and K bands should be 2.3 times higher 
than the values derived from the 2MASS luminosity functions. In Fig. 5, we reproduce the 
luminosity density plots by adding our points. Our K-band luminosity density is about 2 
times higher than those derived from the 2MASS samples(Cole et al. 2001; Kochanek et al. 
2001), much closer but still lower than what the model predicts. The luminosity density of 
our subsample with z < 0.1 is consistent with those from the 2MASS sample. This is not 
surprising since they are derived using similar luminosity functions. Our luminosity density 
for the K < 15 sample is measured at a mean redshift of z ~ 0.1, comparable to SDSS and 
2dFGRS optical samples. However the 2MASS-(SDSS,2dFGRS) paired catalogs are much 
shallower, z ~ 0.05, thus it seems that the most likely explanation is that 2MASS is sampling 
a local underdensity in the galaxy distribution. However it is now clear that at a redshift 
~ 0.15 there is no broad mismatch in optical and IR luminosity densities. We also notice 
that the model, the average spiral galaxy SED (Schmitt et al. 1997) plus a constant tail 
(/„ oc z/ -2 ) at short wavelength (Dwek et al. 1998; Wright 2001), does not fit the luminosity 
densities well, in detail, either. 
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5. CONCLUSIONS 

We present a large and highly complete large K-selected redshift sample down to K < 
15, and use it to derive the K-band luminosity function using both non-parametric and STY 
methods. By comparing with previous K-band luminosity functions, our i^-band luminosity 
function has a significantly brighter characteristic luminosity and steeper slope. The slope 
is the same as in the optical determinations. The K-band luminosity density measured from 
our sample is 2 times higher than those measured from the 2MASS redshift sample, the 
largest local fCselected redshift sample. We argue that this deeper survey, of comparable 
depth to the optical and SDSS surveys, is more strictly comparable in luminosity density 
and that we have in fact resolved the discrepancy of the 2MASS survey. We are also able 
to re-produce the the 2MASS fCband luminosity function using similar sunsample with 
z < 0.1. 

We also derive the if-band luminosity functions for both early and later type galaxies. 
The early type galaxies have a bright M* and a flat slop a ~ — 1 , while the later type galaxies 
have a faint M* and a steep slope a = —1.4. A steep slope for later type galaxies is also 
found in the current large optical redshift surveys. 

This research has been supported by University of Hawaii, AAO, Max-Planck Institute 
for Astronomy, and the Smithsonian Astrophysical Observatory. J.-S. Huang thanks both 
Max-Planck Institute for Astronomy and AAO for their financial support during his visiting 
AAO. 
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Fig. 1. — The redshift distribution for the K < 15 sample of totaling 1056 galaxies. The 
median redshift is 0.138, and the maximum redshift is 0.57. 
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Fig. 2. — The K-hand luminosity function derived from our K < 15 sample with the best 
fit Schechter function. We also plot the la contour at lower right corner. 
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Fig. 3. — The i^-band luminosity functions for the early and later type galaxies, which are 
different from each other: the LF of the later type galaxies has a much steeper slope end 
than that of the early type galaxies; the LF for the early type galaxies has a bright M*. 
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Fig. 4. — Our luminosity function is plotted against two 2MASS X-band luminosity functions 
(Cole et al. 2001; Kochanek et al. 2001), and those of Gardner et al. (1997) and Szokoly et 
al. (1998). We also plot our luminosity function for z < 0.1 which is consistent with those 
from the 2MASS sample. 
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Fig. 5. — The optical and near-infrared luminosity densities derived with Q m = 0.3 and 
= 0.7 cosmology. The optical luminosity densities are measured in u*,g*, r*, i*, and 
z* by the SDSS team(Blanton et al. 2001), and in Bj band from the 2dFGRS (Folkes et 
al. 1999). The J and i^-band luminosity densities are derived from the 2MASS luminosity 
functions (Cole et al. 2001; Kochanek et al. 2001), and from our data. The luminosity density 
derived using our subsample at the low redshift bin (blue filled dot) is consistent with the 
2MASS luminosity density. An average spiral galaxy SED (solid line and the one with an 
flat tail at short wavelength (f u oc is~ 2 , Dashed line) as a background model suggested by 
Dwek et al. (1998) and Wright (2001) are also plotted. 
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Table 1 
X-band Luminosity Function 
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= O.O(corr) 


-23.57±0.08 


-1.33±0.09 


0.017±0.002 


4.82±0.44 


H m 


= 1.0 


n A 


= 0.0(uncorr) 


-23.36±0.08 


-1.16±0.09 


0.020pm0.002 


3.88±0.41 


Q m 


= 1.0 


n A 


= O.O(corr) 


-23.41±0.08 


-1.25±0.09 


0.020±0.002 


4.44±0.40 



1 Luminosity density in units of 10 27 erg/ s/ H zj Mpc 3 

2 LF derived without correction for incompleteness 

3 LF derived with correction for incompleteness 
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Table 2 

Previous and Current K -band Luminosity Functions 



Sample 


M» 


QL 


<t>* 2 


if 


number 




Mobasher et al. 1993 


-23.37±0.30 


— 1 .00 i 0.3 


1.12±0.16 


1.96±0.62 


95 




Glazebrook et al. 1995 


-23.14±0.23 


— 1 .04 i 0.3 


2.22±0.53 


3.19±1.07 


98 


K < 17.3 


Gardner et al. 1997 


-23.30±0.17 


-1.00 ±0.24 


1.44±0.20 


2.36±0.48 


532 


K < 15.0 


Szokoly et al. 1998 


-23.80±0.30 


-1.30 ±0.20 


0.86±0.29 


2.90±1.27 


110 


K < 16.5 


Loveday 2000 


-23.58±0.42 


-1.16±0.19 


1.20±0.08 


2.86±0.94 


345 


b 3 < 17.15 


Kochanek et al. 2001 


-23.43±0.05 


-1.09 ±0.06 


1.16±0.10 


2.27±0.21 


3878 


K < 11.25 


Cole et al. 2001 


-23.36±0.02 


-0.93 ±0.04 


1.16±0.17 


1.94±0.29 


5683 


K < 13.2 


This paper 


-23.70±0.08 


-1.39 ±0.09 


1.30±0.20 


4.79±0.16 


1056 


K < 15.0 



1 LF derived with Q m = 0.3 and Q A = 0.7. 

2 4>* in units °f W~ 2 h 3 Mpc~ 3 . 

3 Luminosity density in units of 10 27 erg/ s/ H z/ 'Mpc s 
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Table 3 

E"-band Luminosity Functions for early and later type galaxies 





Early 


Type 


Later 


Type 


Sample 


M» 


a 


M, 


a 


Kochanek et al. 2001 


-23.53±0.06 


-0.92 ± 0.10 


-22.98±0.06 


-0.87 ±0.09 


This paper 


-23.56±0.26 


-1.03 ± 0.31 


-23.28±0.28 


-1.42 ± 0.31 



